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ABSTRACT 

Analysis of nearly 900 bathythermograms taken at half- 
hourly intervals at Ocean Weather Station ECHO between 
2 and 21 September 1959 reveals that convective mixing, 
the most easily identified physical process involved in the 
daily heating and cooling cycle, extended to a depth of 60 
feet near sunrise. Measurement of radiational heating was 
complicated by apparent internal wave action on a weak 
thermocline at 60 feet. Estimation of the coefficient of 
thermal conductivity from the data yielded a value much 
lower than expected because of simultaneous turbulent 
heat transfer. After elimination of known sources of tem- 
perature variations (instrumental or human errors), advec- 
tive heat exchange was revealed to be the major source of 
remaining variations. The energy involved in advection was 
much greater than that involved in the vertical processes 
of heating, cooling, and convective mixing. 


FOREWORD 


Short-term prediction of environmental factors: is an essential 
part of the AntiSubmarine Warfare: Pavironmential Prediction System 
(ASWEPS). The physical processes affecting oceanic environments must 
be understood before prediction techniques can be formulated. Recog- 
nizing that observations of diurnal variations in the open ocean are 
scarce, the Oceanographic Office has conducted a series of observation 


programs covering a Limited area in cooperation with the U. S. Coast 


Guard. 
This report describes one aspect of the observational results: 


an analysis of the diurnal processes in the mid-Atlantic during sum- 


mer. Therefore it serves as a forerunner to a complete understanding 


of the processes which, in general, govern short-term environmental 


changes. 
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INTRODUCTION 


The ASWEPS research program being conducted at the Hydrographic Of- 
fice* includes a systematic investigation of the physical processes in- 
volved in the daily heating cycle of the upper layers of the ocean. This 
investigation is being conducted in order to determine the effect of 
these processes on the thermal structure and to devise an analytical 
short-range prediction technique. The approach follows the method out- 
lined by Schule (1952). A thermal structure prediction program based on 
meteorological parameters has been tested by this office (Nix 1962). 
This report summarizes the results of an attempt to find the numerical 
values of some vertical physical processes in the water under conditions 
of minimal horizontal variations. 


Many studies of diurnal temperature variations in lakes have been 
made (for example, Geiger, 1950); fewer studies of diurnal changes in the 
oceans have been made (for example, Sverdrup et al., 1942 and Defant, 
1961). Perhaps the most extensive account is given by La Fond (1954) who 
considered inshore as well as deep-sea locations. Most of the oceanic 
studies were made at anchor stations of short duration. Because they are 
limited to shallow water, a comparatively large number of studies have 
been excluded from this discussion. 


This report analyzes the results of an intensive survey conducted by 
personnel of the Hydrographic Office in September 1959 at Ocean Weather 
Station ECHO (35N, 48W). Half-hourly deep bathythermograph (BT) casts 
and several short series of shallow casts at 10-minute intervals were made. 
In addition, 93 Nansen casts were made in connection with a study of the 
seasonal halocline. BT data were carefully controlled for accuracy by 
reading slides immediately upon completion of an observation; instruments 
were replaced when gross errors appeared in the data. Thirteen instru- 
ments were used to obtain nearly 900 BT’s between 2 and 21 September 1959. 
The mean diurnal temperature changes for a mid-latitude mid-oceanic loca- 
tion in late summer can be obtained from these observations. 


METHOD OF ANALYSIS 


All BT records were averaged for each hour of the day. Thus, the 
average for each hour represents between 27 and 49 observations taken dur- 
ing the period of the survey. The purpose of averaging is twofold; namely, 
(1) to provide hourly values of meteorological and astronomical variables 
for the computer prediction program, and (2) to remove or minimize hori- 
zontal factors, mainly advection and the horizontal component of diffusion, 
which are not used in the computer program. In addition, some of the ver- 
tical processes average out in desirable fashion; that is, the average 
thermal structure represents averages of cloudiness, wind (evaporation), 
sea condition (and therefore average turbulent mixing), and local cooling 
due to shower activity. Since the observations were taken over a period 
of 20 days, during which times of flood and ebb tides moved forward about 
17 hours, tidal effect has almost been removed from the data. 


* Redesignated U. S. Naval Oceanographic Office 10 July 1962 


SURFACE CHARACTERISTICS 


In order to determine whether seasonal influences were important, 
the data were divided into two time periods, covering approximately 10 
days each. Mean hourly sea surface temperatures are plotted in Figure 
1, which provides a graphic comparison of the mean temperature varia-= 
tions of the periods 2 to 10 September and 11 to 21 September. The 
curves are similar in shape and have almost identical mean values. Al- 
though slightly more heating occurs in the afternoon hours during the 
earlier period, data for the two periods were combined with the subsur- 
face temperatures for the same period. The curves are nearly identical 
in value at night. The mean daily range in sea surface temperature was 
about 1.1° F. 


SUBSURFACE TEMPERATURES 


Subsurface temperatures were Listed according to the standard SERC 
deck coding procedure used by this Office. This procedure records tem- 
peratures at 20-foot intervals between the surface and 360 feet and at 
hO-foot intervals between 360 and 880 feet. ‘The influence of the daily 
temperature cycle was not expected to extend below a depth of 100 feet; 
therefore the mean hourly temperatures were computed only for the sur- 
face and for every 20-foot level to this depth. Table 1 lists the 
hourly means and standard deviations. The means have been plotted in 
Figure 2, which also shows times of local sunrise, noon, and sunset. 


The temperature apparently fluctuates rather regularly and con- 
siderably at all levels; separation of the heating and cooling effects 
from those due to other influences such as internal waves becomes a 
problem. This problem was solved by analyzing the thermal structure 
in layers. Major daily thermal variations result from incoming solar 
radiation during daylight hours and from loss of heat during night hours. 
The latter of these effects may result in convection if density insta- 
bility is produced. 


As an aid to analysis, the temperature differences in five layers 
have been plotted (Figure 3). The resulting curves indicate that con- 
vection (zero gradient) began at 20 feet at 02002, LO feet at 05002, 
60 feet at 1000Z, and that it did not extend to 80 feet. 


Figure 4, a further exposition of the layer analysis, shows the 
mean hourly temperature differences between the surface and 60 feet as 
compared to those in the layer between 60 and 100 feet. Diurnal changes 
are immediately apparent from the regularity of the curve for the upper 
layer and lack of regularity of the curve for the lower layer. There- 
fore, it can be concluded that diurnal effects were confined to depths 
between the surface and 60 feet at OWS ECHO in September 1959. 
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Mean Hourly Temperatures and Standard Deviations (° F) 


Surface 


20 ft 


ho ft 


60 ft 


80 ft 


77-33 
6 


100 ft 


No. of obs. 


KL 


896 


YySEGWaldas Ie 


L3SNNS 
1v907 


2 ‘SAYUNLVYSdWAL ATHNOH NVAW 2 3yNdIS 


oe Ee 


Eas ire 
PTT | avoon | ET | voor [| TTT 


(4 0) SYNIVYSdIW3L 


S13A31 39VSYNSENS GNV 39VsYNS 
N3SML38 SSONSYSSSIG SYNIVYUSUNSL AITYNOH NVAW € 3YNDIS 


Z) JWI 


C-0- 
si HIS: sna aca 


Nig TN. EEL AE 
kag st at Fea Dc 
PWN Vets e Saas eee 
JAB OEn Cee Mee Vinee aes 
I e972 ere ea a a 
[Se Ga les eal UE ee ee SU 
HHA ee EEE 


£0 


vO 


To) 
fo) 
do) JIN3IY3SII0 SERENE. 


ZO000! lv 13 09 
Z00SO lv 14 0¥ 
20020 1V 14 O2 
OY3Z SI LN3IQVYS SHL N3HM SYNIDIO NOILISANOD 


EO 


(4 5 a aN YG ayereticrint: 


tk 
H tad 
wi} 
, wt 
(e} 
wo 
te 
} << 
a} 
Ww | 
Qc 
Oc § 
DI} 
oO; 
OF; 
oO; 
= 
jr 
Ee] 
oO 
} Ud 
1 >| 
} x 
oO 
(ey 


TIME (Z) 
FIGURE 4 HOURLY TEMPERATURE DIFFERENCES IN TWO LAYERS 


The above conclusion helps to explain hourly temperature fluctua- 
tions at 60 feet. The two regular wave cycles between sunset and sunrise 
in Figure 5 are apparently caused by internal waves. Two fluctuations 
during the daylight hours are partly masked by diurnal heating. If the 
curve for the 60- to 100-foot layer in Figure is correct, heat changes 
at 60 feet were directly reflected at 100 feet. ‘The maximum temperature 
at 17002 at 60 feet (Figure 5) was about 0.2° F higher than the secondary 
maximum at OOO0OZ. This difference may be interpreted as the mean diurnal 
temperature range at 60 feet. Rough correction for this value in Figure 
5 fails to eliminate the two fluctuations during the daylight hours. 
Thus, four internal waves are apparently revealed by the hourly analysis. 
The analysis automatically eliminates waves with frequencies of less than 
2 hours. Since tidal effects have been largely eliminated by averaging 
the data over a period of 20 days, these internal waves must be station-~ 
ary and are probably caused by some other phenomena. One theory of a 
possible cause is isostasy; that is, the waves may be caused by convec= 
tion during which the water particles seek equilibrium after vertical 
transfer from the surface. Inasmuch as the convective process required 
3 hours to progress from 20 to 40 feet and 5 hours to progress from ho 
to 60 feet, this theory appears to be an unsatisfactory explanation of 
- the internal waves. 


Figures 6, 7, and 8 show the results of the convective process. 
Figure 6 illustrates depth of convection and time of maximum heating. 
Figure 7 shows mean BT's for each hour during the night; Figure 8 shows 
the same for daylight hours. 


SENSIBLE HEAT CHANGES 


Temperature fluctuations in the upper 60 feet are analogous to, but 
not identical to, changes in sensible heat content, which depend on the 
density changes. Salinity observations were obtained from the 93 Nansen 
casts. Possible diurnal changes in salinity could not be determined be=- 
cause these observations were not evenly distributed throughout the day. 
The majority were taken at 6-hourly intervals (03002, 0900Z, 15002, and 
21002). Mean salinity values in Figure 9 show little diurnal vertical 
variation in the upper 20 meters; whereas the diurnal range of salinity 
was considerable at 40 meters. This range is probably due to fluctua- 
tions of the seasonal halocline. A mean salinity value of 36.25 °/o0 
was used with standard formulas for converting temperature into sensi- 
ble heat. 


Figure 10 shows hourly changes of total sensible heat content (Q) 
in the three 20-foot layers of the ocean. The curves closely resemble 
those of Figure 2, owing to use of the constant salinity value in the 
formulas. In order to outline the difference in heat content between 
layers more clearly, the differences between the first and second layers 
and between the second and third layers -have been plotted in Figure 11. 
The relationship between heat content changes caused by solar heating 
and convection is also shown. The fact that increases in heat content 
due to solar radiation were practically simultaneous in all layers should 
be noted in this Figure. This simultaneity is due to high transparency 
of the water. Cooling due to heat loss at the surface progressed from 
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FIGURE 5 MEAN DIURNAL TEMPERATURE FLUCTUATIONS AT 60 FEET 
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FIGURE 6 PHYSICAL PROCESSES OF THE DIURNAL HEATING CYCLE 
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the first layer to the third layer. In summation, the maximum heat 
differential between the first and second layers and between the 
second and third layers did not occur until 18002 and 200072 respece- 
tively, or nearly at sunset, although the maximum surface tempera- 
ture was reached at 1700Z or 2.3 hours after local noon. 


Figure 12 continues the study of differential heating of layers 
and shows that the differences of heat content between the third and 
fourth layers and between the fourth and fifth layers do not show 
much regularity. This indicates that the heating process is confined 
. mainly to the upper 40 feet, although convection extends to 60 feet 
if differential heating alone is considered. In contrast, reversed 
Logic leads to the conclusion that appreciable solar heating actually 
extends well below 100 feet as shown by the peaks in Figure 2. Howe 
ever, the amount of radiation absorbed in each layer is more constant 
than would be assumed from normal extinction coefficients (Sverdrup 
et al., 1942). As explained above, it must be remembered that Figure 
2 is complicated by apparent internal wave action below a depth of 60 
feet. 


COEFFICIENT OF THERMAL CONDUCTIVITY 


An inherent property of substances is equalization of temperature 
within the substance (assuming absence of outside influence). This 
property is thermal conductivity and occurs when thermal gradients ex- 
ist in the substance. The coefficient of thermal conductivity for sea 
water is slightly smaller than that for fresh water, according to the 
Smithsonian Tables. The BT data described above may be used to esti- 
mate this coefficient. 


Since the coefficient of thermal conductivity is given in units 
of calories per centimeter per second per degree centigrade, it is evi- 
dent that the dimensions of these units are equivalent to the rate of 
change of the temperature gradient. If this rate of change can be ap- 
proximated, a value of the coefficient can be determined. In terms of 
finite differences, the temperature difference between layers is pro- 
portional to the temperature gradient. The time change of the differ- 
ence is the first differential, and the rate of change is the coeffi- 
cient of thermal conductivity. 


Sampling errors are minimized by use of fitted curves. Therefore 
a Fourier analysis (for the second harmonic) of the differences plotted 
in Figure 3 was performed by digital computer. The fitted curves are 
shown in Figure 13. Gradients within the layers are shown in Figure lh. 
Table 2 lists changes in temperature differences, from which the coeffi- 
cient may be estimated. In Table 2, phase changes from positive to neg- 
ative or differences of maximal differentials are about one hour between 
the 20- to 40-foot layer and the 40- to 60-foot layer. ‘The O- to 20- 
foot layer must be omitted because of wave action. As shown in Figure 
14, temperature gradients in the upper two layers (0-20 and 20-40 feet) 
are in phase; however, wind mixing reduces the gradients in the 0O- to 
20-foot layer so that eddy conductivity and thermal conductivity cannot 
be distinguished. 
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TABLE 2 


Change in Temperature Differentials from Fitted Curves (° F) 


Hour (Z) 0-20 ft 20-40 ft 40-60 ft 60-80 ft 
1200-1300 2022 +026 2008 2005 
1300-1400 2027 0032 0012 2006 
1400-1500 2028 2035 018 00} 
1500-1600 2026 2033 2022 OO} 
1600-1700 02% 2028 - .022 200k 
1700-1800 e012 2018 2020 e001 
1800-1900 2002 2006 2014 ~000 
1900-2000 ~ 008 -.007 ie? 2005 = 2002 
2000-2100 - O17 -.017 - 00} = 006 
2100-2200 -.022 = .022 Sorat -.007 
2200-2300 ~ 024 -.028 -.021 -.010 
2300-0000 -.023 = 2027 -.023 -.002 


If the rate of phase difference is estimated to be 20 feet per 
hour, the amount of change may be estimated. The average absolute tem- 
perature differential change within the 20- to 40-foot layer in any 
given hour and within the 40- to 60-foot layer an hour later is .0073° F. 
Using a specific heat of .93 cal/gm/° C, the thermal conductivity co- 
efficient is .00065 cal/cm/sec/° C. The figure of .00065 is considerably 
lower than that given in the Smithsonian Tables (00143). Thus some 
fifty percent of thermal conductivity appears to be counterbalanced by 
turbulent mixing due to wave action. It should be noted that Table 2 
includes only daylight hours during which convective mixing was absent, 
so that turbulent or mechanical mixing constitutes the only process for 
vertical heat transfer that is unaccounted for by this analysis. 


VARIABILITY 


In the discussion above, conclusions have been confined to the ex- 
planation of mean temperature and mean salinity changes in the upper 
layer of the ocean. However, as with all observations, there is much 
variability in the temperature readings from which the mean values are 
computed. The major sources of variation may be summarized as follows: 
(1) Instrumental error; that is, errors caused by faulty temperature and 
pressure elements in a bathythermograph.e Such errors may be assumed 
constant for temperature and to increase slightly with depth as pressure 
increases. (2) Reading error, defined as the average difference between 
readings of the same BT trace by two different persons. (3) Interdiurnal 
variability; that is, day-to-day difference in temperature owing to either 
short- or long-period fluctuations such as tidal periods and seasonal ef- 
fects, plus variability introduced by grouping the observations by hours. 
(4) Advection or horizontal movements of water of differing temperatures. 
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Although the actual values of these variations are not known 

with certainty, tentative values may be given. A study of scale error 
in BI's made by this Office in 1959 revealed an error of approximately 
0.2° Fe. The reading error is estimated at about 0.19 F. The inter- 
diurnal variation is small, probably no more than 0.29 F, but is con- 
centrated near the surface because of wave action. Thus, after allow- 
ing for these sources of variation, Table 3 shows that the major source 
of variation must be presumed to be advection. 


TABLE 3 


Hourly Standard Deviations of Temperature After 
Removal of Variations Other than Advection (° F) 


Hour (Z) Surface 208 ho ft 60 ft 80 ft 100 ft 
0000 059 269 °72 83 °19 91 
0100 5 Sif oS 261 68 i 82 
0200 257 058 °6L 63 68 073 
0300 059 061 263 -66 ar) oT4 
0400 259 059 59 Alon 265 + 669 
0500 065 065 Ok 067 -69 82 
0600 62 261 063 63 073 88 
0700 oA 5 346 046 alli 251 o5T 
0800 558) oS 058 (Sit °61 63 
0900 058 057 256 555 052 62 
1000 66 268 68 269 Bri 286 
1100 69 268 69 ie. 82 98 
1200 268 269 Se aD 22 092 
1300 69 073 Th °78 sO6L iam 9k 
1400 053 °60 ~60 065 72 81 
1500 258 059 058 263 RAS ok 
1600 Soil oh6 oS oil ait °72 
1700 058 sgh “55 60 oS peal 
1800 aval, 065 66 68 oh e81 
1900 259 OBS 62 ail 076 86 
2000 067 065 66 078 oOL 283 
2100 015 °70 -67 076 82 088 
2200 Onl nt °72 73 °76 283 
2300 Her it rcs °76 73 083 


A standard deviation of 0.1° F is equivalent to 130 gm cal/cem in 
a layer 20 feet thick. In Figure 3 it can be seen that the temperature 
gradient for the top 20 feet decreased from 0.1° to 0° F between 2300Z 
and 0200Z, or three hours, corresponding to a heat change of 22 gm 
cal/em@/hr. The gradient in the upper 40 feet also changed from 0.1° 
to 0° F between 0100Z and 0500Z, or four hours, a rate of change of 32 
gm cal/cm-/hr. Similarly the gradient in the top 60 feet was reduced 
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to 0° F between 0300Z and 10002, or 7 hours, a rate of 28 gm cal/cm@/hr. 
These calculations are summarized in Table 4. These values are very 
small compared to the advective variation which amounts to at least 

590 gm cal/cm@/hr for the smallest standard deviation in Table 3. How- 
ever, the vertical changes are much more regular and also much better 
organized owing to the physical processes involved. 


TABLE 4 


Energy Involved in Convective Mixing 


Hour (Z) Depth of Convection Energy 
(em cal/cm?/hr ) 

2300-0000 2{0) he 22 
0000-0100 20 ft 22 
0100-0200 ho ft 54 (22 +32) 
0200-0300 ho ft 32 
0300-0400 60 ft 60 (32 +28) 
0400-0500 60 ft 28 
0500-0600 60 ft 7 28 
0600-0700 60 ft 28 
0700-0800 60 ft 28 
0800-0900 60 ft 28 
0900-1000 60 ft 28 

SUMMARY 


This study was undertaken for the purpose of estimating the nu- 
merical values of various physical processes involved in the diurnal 
eycle of heating and cooling in the upper layer of the ocean. Magni- 
tude of the vertical processes, convection,and surface heating, hagvq 
been estimated with some degree of precision. However, the energy 
contributed by the major horizontal process of advection is much 
greater than that contributed by the vertical processes. In order to 
incorporate advection into the deductive and predictive procedures 
which are under development for the ASWEPS program, either a greater 
number of accurate, closely-spaced synoptic observations must be ob- 
tained, or thermal structure prediction must be confined to long time 
periods in which advection can be minimized. Otherwise, direct means 
of measuring advection must be devised. 
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